On the size of linear superpositions in molecular nanomagnets 
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Molecular nanomagnets are relatively complex spin systems that exhibit quantum mechanical 
behavior at low temperatures. Exploiting quantum-information theoretic measures we quantify 
the size of linear superpositions that can be generated within the ground multiplet of high-spin 
nanomagnets. In particular, we consider the prototypical single-molecule magnets (namely Mni2 
and Feg) and a physically significant class of spin rings (including Mng), all characterized by a 
ferrimagnetic spin ordering in the ground state. These cases are compared with those achievable 
in ferromagnetic spin clusters. We show that the size of these linear superpositions can be further 
enhanced by increasing the asymmetry between the sublattices, and by reducing the competition 
between exchange interactions within the nanomagnets. 

PACS numbers: 03.65.Ta,75.50.Xx 



Molecular nanomagnets provide prototypical examples 
of zero-dimensional spin systems, with physical proper- 
ties that can be widely tuned by chemical synthesis lj. 
Antifcrromagnetic exchange interactions typically result 
in the presence of different forms of quantum entangle- 
ment at low temperatures, such as those between individ- 
ual spins of a single cluster [2] or the collective spins of 
weakly coupled molecules [3J. Quite generally, the inter- 
est in molecular nanomagnets was aroused by the early 
observation of quantum effects. These include quantum 
tunneling of the molecule spin, resulting in a speed-up of 
the magnetization relaxation [?J [5] , and quantum phase 
interference [5J . More recently, microwave-induced quan- 
tum coherences [5] and oscillations HO] were ob- 
served in a wide class of spin clusters. These experi- 
ments demonstrate the possibility of manipulating the 
molecule spin state on time scales shorter than the deco- 
herence time, which represents a preliminary requirement 
for the implementation of quantum-information process- 
ing nn na. 

In all these cases, the relevant quantum coherences 
are those between states with different spin projection, 
belonging to the ground S multiplet of the molecular 
nanomagnct. These linear superpositions can be partic- 
ularly meaningful in high-spin molecules, also referred to 
as single-molecule magnets. Here, the ground S mul- 
tiplet may include classical-like states that are meso- 
scopically different from one another, and whose lin- 
ear superpositions can thus be regarded as mesoscopic 
Schrodinger cat states [T3J. Hereafter we quantitatively 
investigate the size of such linear superpositions, start- 
ing from realistic spin Hamiltonian models of the nano- 
magnets. We specifically refer to the most celebrated 
single-molecule magnets, [Mni20i2(CH 3 COO)i6(H 2 0)4] 
HD and [Fe 8 2 (OH) 12 (tacn) 6 ]Br 8 (H 2 0) 9 (tacn=l,4,7- 
triazcyclononane) |15j , and to a class of spin rings whose 
prototype is represented by the [Mn(hfac) 2 (NITPh)] 6 
molecule (Mug) [15]. In all the three molecules (here- 
after refereed to as Mni 2 , Fes, and Mng) the ground 



state is characterized by a ferrimagnetic ordering, with 
the total spin S resulting from the inequivalence of the 
two sublattices [1] . The size of possible linear superpo- 
sitions generated within the ground multiplet is favored 
by the absence in Mng of competing interactions within 
each sublattice, and by the large unbalance between their 
partial spin sums. In Mni 2 and Fes, the more complex 
pattern of partially competing exchange interactions lim- 
its the classical-like character of the ground states, and 
thus the measure of the linear superpositions that takes 
this aspect into account. 

The problem of quantifying the size of a linear superpo- 
sition has been addressed from different perspectives, in 
order to assess the actual macroscopicity of quantum co- 
herences in infinite physical systems [I7H22] . There, the 
question whether or not a linear superposition is macro- 
scopic is answered by considering the limit of different 
measures as the number N of the microscopic subsys- 
tems tends to infinity. In the present paper, we use two 
of such measures to quantify the size of linear superpo- 
sitions in zero-dimensional spin systems. The first mea- 
sure is based on the quantum Fisher information [19] . 
and quantifies both the non-classicality of the linear su- 
perposition |*) = (|*i) + |* 2 ))/\/2 and the classical- 
like character of its components in terms of the respec- 
tive quantum fluctuations of single-spin operators. As 
showed in following, the advantage of such measure in the 
present context is at least twofold. Firstly, it quantifies 
the distinguishability between the two superposed states 
by means of a physically intuitive quantity, such as the 
spin projection of the constituent ions. Secondly, it can 
be expressed in terms of spin-pair correlation functions, 
that are now experimentally accessible in molecular spin 
clusters [23] • The second measure we consider quantifies 
the size of a linear superposition in terms of the possi- 
bility to distinguish between the two components |^i) 
and l^) by means of local measurements. The size of 
I ty) is thus identified with the number of subsystems into 
which the spin cluster can be partitioned, such that the 






FIG. 1: (color online) (a) Spin cluster corresponding to the Fes molecule. The eight s = 5/2 spins are divided into the two 
sublattices A (light and dark blue) and B (red), defined by their relative orientation in the system ground state. Equivalent spins 
are denoted by the same color. Same line style corresponds to identical values of the exchange coupling between neighboring 
spin pairs, (b) Spin cluster corresponding to the Mni2 molecule, where s t = 2 and sf — 3/2, with the same conventions as 
above, (c) Schematics of the double-well potential induced by the crystal field within the ground S = 10 multiplet of the two 
spin clusters, and corresponding level scheme. 



which-componcnt information is available in each of the 
subsystems [24] • The notion of distinguishability can be 
itself translated in terms of measurement outcomes, but 
in a way that makes it of no experimental relevance in the 
present context. In fact, the projective measurement of 
the individual spins (or of arbitrary subsets of them) that 
would be required to estimate the distinguishability are 
presently unfeasible in molecular nanomagnets. On the 
other hand, this measure has the advantage of being inde- 
pendent on the specific class of one-body operators that 
we use with the Fisher information, and is thus suitable 
for comparing the size of linear superpositions obtained 
in completely different systems. 

The paper is organized as follows. In Section I we in- 
troduce the measures that are used to quantify the size 
of the linear superpositions. In Section II we present the 
results obtained for the Mni2 and the Fes spin clusters. 
In Section III we discuss the case of spin rings, includ- 
ing the Mng molecule. Section IV and V are devoted to 
a comparison with other molecular spin clusters and to 
the conclusions, respectively. Finally, in Appendix [A] we 
prove the optimality of the staggered magnetization as a 
one-spin operator for the Fisher information. 



I. SIZE OF LINEAR SUPERPOSITIONS IN 
SPIN CLUSTERS 

We consider linear superpositions of the form 
1 



I*', =-^(|*i) + |* a )), 



(1) 



where |^i) and l^) are two eigenstates of the spin 
Hamiltonian Hq of the molecular nanomagnet, that com- 
mutes with the total-spin projection S z . More specifi- 
cally, we focus on the case where the dominant terms of 
H are represented by the exchange couplings between 
the N spins of the cluster, and the states |\&i) and j^) 
belong to the ground iS-multiplet, with different values 
of S z (Mi and M 2 , respectively). 



A. Quantum Fisher information 

The first criterion we use for quantifying the size of the 
linear superposition is based on the quantum Fisher in- 
formation [T5]. For pure states, such quantity coincides, 
up to a multiplicative constant, with the variance of the 
relevant operator X: 



F*(X)=4[(X*)-(X) 2 ]. 



(2) 



The one-particle operator X that maximizes J-V should 
coincide with the one that is best suited to distinguish 
between the two components. The operator X is thus 
expected to have largely different expectation values for 
l^i) and \^2), and small fluctuations within each com- 
ponent. In the class of systems considered hereafter, this 
criterion is fully met by the staggered magnetization op- 
erator St: 
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where A and B are two sublattices into which the spin 
cluster is partitioned (Na + Nb — N). For each sys- 
tem considered hereafter, each sublattice is defined so 
as to include all the spins Sj with the same sign of the 
expectation value (^k\sz,i\^k}> for any given component 
k = 1,2 of the linear superposition. The staggered mag- 
netization S* maximizes the Fisher information for the 
most relevant linear superpositions of the form Eq. [I] 
(namely, those between components with large values of 

| Mi — Mi | ) , within the class of operators X = J2i=i **i ' s » 
(with nf = 1) (see Appendix [A[) . Such linear combina- 
tions of single-spin terms have an intuitive physical mean- 
ing and represent the most significant operators from an 
experimental viewpoint, for their variance can be decom- 
posed into experimentally accessible quantities, such as 
spin-pair correlation functions |23j . Quadratic functions 
of Si ia (with a — x,y, z), are instead unsuitable for char- 
acterizing the above linear superpositions, for they give 
identical expectation values for pairs of states that are 
the time-reversal of one another. Cubic (or higher) pow- 
ers of the single-spin operators don't correspond to any 




FIG. 2: (color online) (a) Effective size Dfi^) of the linear 
superpositions j*) formed by the ground states of the Mni2 
molecule with total spin projections Mi and M2. The maxi- 
mum possible value for the spin cluster is J]V Si — 22. In the 
cases Mi = M2 we set Dfi = 0. (b) Size of the same lin- 
ear superpositions based on the relative Fisher information, 
Drfi(V). 



observable quantity in molecular nanomagnets, and even 
less do their fluctuations pQ. 

For a given spin cluster, the maximum of the Fisher in- 
formation corresponding to the staggered magnetization 
is given by: 
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where the norm of the spin operator is here defined as its 
largest eigenvalue (in modulus). This value corresponds 
to the state |tf rooa! ) = (|*™ aa: ) + \^™ ax ))/^/2, where the 
components are given by 
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(5) 



As discussed below, these states can also be regarded as 
the most classical-like states of the spin cluster, and have 
vanishing fluctuations of the staggered magnetization op- 
erator. 

The quantum Fisher information can be used to define 
the effective size of the linear superposition, along the 
lines of Ref. [T5]: 
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(6) 



Such quantity can take any value in the interval 

[0, J2y=A b J2i=i s f]- The maximum size of a linear su- 
perposition that can in principle be generated in a given 
spin cluster thus depends not only on the number N of 
the constituent spins, but also on their values Si- This fol- 
lows from the use of spin operators - rather than generic 
norm-one operators - in the quantum Fisher information. 
It also assigns a larger size to the linear superposition 
of the states |\£j£?f 2 ), with respect to, e.g., one where 
the components are given by the product of single-spin 
states \mf k ) (x — -^ B) with non-maximal spin projec- 
tions (|m* fe | 7^ sf). We note that such discriminations 
also allows to link the present measure Dpj^) to the 
ones that quantify the size of a linear superposition in 



terms of its decoherence rate [20]. In fact, the fragility 
of a linear superposition between two eigenstates of a 
molecular spin cluster was shown to depend on their dis- 
tinguishability in terms of single-spin expectation values 
[2 5) , and is thus an increasing function of the quantities 
|<*i|SSI*i) - (* a |SJ|* 2 )|. 

The absolute value of the Fisher information charac- 
terizes the linear superposition in terms of the quantum 
fluctuations of the operator S 1 *. However, it doesn't dis- 
criminate between the contribution to such fluctuations 
coming from the components \^i) and ^2); an d the one 
that results from their linear superposition. In the spirit 
of the Schrodinger gedanken experiment |13j , the super- 
posed states should in fact approach a classical charac- 
ter, and therefore present the smallest amount of quan- 
tum fluctuations allowed by the uncertainty relations. In 
this perspective, a suitable figure of merit for the linear 
superposition is represented by the relative Fisher infor- 
mation, where the variance of the relevant operator in 
|$) is normalized to that in the components J19) : 



D 
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J* (S*) 
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(7) 



For a spin cluster, the nonclassicality of a state can in 
general be traced back to the presence of entanglement 
between the spins and of quantum fluctuations in the 
single-spin operators. A classical-like \^k) can thus be 
identified with the product of single-spin coherent states 
|0), defined by the relation |(</>|s|0)| = s [25]. In the 
present case, the only relevant coherent states are repre- 
sented by \m — ±s), due to symmetry reasons. In fact, 
being the two components of Eq. [T] eigenstates of S z , the 
single-spin reduced density matrix is given by 



o {1) 

Pi,k 



tr j¥i {|tt*)<* fc |} 






p fc (m i )|m i )(rn l 



(8) 



and the in-plane components of the spin have vanish- 
ing expectation values: (^k\s i}X \^k) = (*fc|si,y|*fc) = 
0. The classical-like character of the states \^k) is 
thus quantified globally by the fluctuations in the stag- 
gered magnetization (J r ^, fc (S'*)), and, at the single- 
spin level, by the von Neumann entropies S(.Si,k) — 
~ Y^ mi Pk(mi) ln[pfc(m;)] or by the populations Pk(rrii = 
±Sj). In order to maximize the above sizes of the lin- 
ear superposition |\&), the components should approach 
the form \^ ax ) (Eq. [If, characterized by Fv^S*) = 
S(si, k) = 0, and pfc(ra, = ±s») = 1. 

We finally note that general symmetry arguments al- 
low to simplify the calculation of the above quantities. 
In the relevant case where the two components are the 
time reversal of one another (Mi = — M2), the following 
relations apply: 

(*i|STI*i) = -(*2|s;|* 2 ), ms*M = (9) 

(* fe |(5:) 2 i* fc ) = (*i(^) 2 i*>- (10) 

The expression of the relative Fisher information thus 
simplifies to: 
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FIG. 3: (color online) Probability P ; of discriminating be- 
tween the two components l^i) and |^2} (Mi = —M2 = M) 
of the Mni2 spin cluster. The probabilities refer to subensem- 
bles formed by q spins: blue, green and red squares correspond 
to the spins s 1 , s 2 and sf , respectively (q = 1); the purple 
squares correspond to the pair formed by the spins s 1 and sf 
(q — 2). Inset: Values of the von Neumann entropy S(sf,k) 
for the three inequivalent spins corresponding to each compo- 
nent \^k)- Same convention as above for the colors. 



More generally, the values of the Fisher information cor- 
responding to different linear superpositions of the above 
type (Eq. [I]) are not independent from one another. In 
fact, the operators S* and (S*) 2 can be written as com- 
ponents of irreducible tensor operators, whose expecta- 
tion values within a given S'-multiplet are related by the 
Wigner-Eckart theorem [27]. In particular, S* — S^, 
is the component q = of a tensor with rank p = 1, 



whereas (S*) 2 



s^/Vs- 



Sq 'y/2/y/Z can be writ- 



ten as the sum of a scalar and of the 0-th component 
of a tensor with p = 2. As far as S* is concerned, the 
Wigner-Eckart theorem thus implies: 



(aSM\S* z \aSM) = (aS\\S ( 



Ml 



aS)A M0M , 



(12) 



where A^ j2 ^ m are the Clebsch-Gordan coefficients and 
a denotes the ground S'-multiplet of the spin Hamilto- 
nian. The expectation values of (S*) 2 obey instead the 
equation: 



(aSM\{S* z Y\aSM) 



(aS||S< 2) ||o 



' <2S V2/V3. (13) 



S)Amom 



In the following, it will thus be sufficient to compute 
the values of S* and (S*) 2 for one and two different val- 
ues of M, respectively, in order to derive all the others 
from the reduced matrix elements. The same applies 
to the single- and two-spin expectation values discussed 
throughout the paper. 



B. Distinguishability by local measurements 

The effective size of the linear superposition |W) can 
also be quantified in terms of the maximal number of 



subsystems that allow to distinguish between |>I>i) and 
1^2) with probability higher than 1 — 6 [24]. The max- 
imum probability of successfully discriminating the two 
states by performing an optimal measurement on q spins 
of the cluster is given by: 



P 



(9) _ 



1 1, 



(«) 



pSI 



(14) 



,(«) 



where p\ I is the reduced density matrix of the ^-th sub- 
system formed by q spins, obtained from the state |\tfc); 
||X||i := 5^,- 1 Ai| (with A,; the eigenvalues of X) is the 
trace norm. For the TV-spin clusters considered hereafter, 
the effective dimension based on the distinguishability of 
the two components by means of local measurements is 
given by: 



D L m(^,S) = maxn, 



(15) 



where n is the number of subsystems Ai into which the 
spin cluster is partitioned, and the maximization is per- 
formed over all the partitions such that P; > 1 — 6 for all 
I = 1, . . . , n, with 1J™ =1 Ai = S and S the TV-spin system. 
The size quantification based on the local distinguisha- 
bility has the advantage, with respect to the previously 
defined measure, of being independent on the specific 
class of operators X. It thus allows to compare the 
sizes of linear superpositions obtained in different sys- 
tems. On the other hand, the connection to quantities 
that are experimentally accessible in molecular nanomag- 
nets is not straightforward. The maximum probability 
P{ is obtained by a projective measurement in the basis 
that diagonalizes pfl — p\ g (see Ref. [24] and references 
therein). In the case where the subsystem Ai is repre- 
sented by a single spin Sj, given the form of the reduced 
density matrix (Eq. ^), the optimum measurement is 
that of s z ,. In the case where A\ includes more than 
one spin, the expression of the relevant observable can- 
not be determined a priori, and might be a non-trivial 
function of the spin operators. In any case, projective 
measurements of single spins, or of specific subsystems 
of the spin cluster, are presently unavailable in molecu- 
lar magnetism. 



II. SINGLE-MOLECULE MAGNETS 

The prototypical high-spin molecules are probably rep- 
resented by Mni2 and Fes [I]- Hereafter we compute the 
size of the linear superpositions that can be obtained by 
linearly combining cigenstates of their spin Hamiltonian 
belonging to the ground S multiplct. 

The dominant part Hexc of the spin Hamiltonian of 
both these molecules corresponds to the isotropic ex- 
change interaction between neighboring spins. There- 
fore, the total spin S can be regarded as a good quantum 
number. Smaller terms, such as those related to crystal 
field, can be treated at a perturbative level and projected 
within the ground S'-multiplet 1 . In particular, the crys- 
tal field anisotropy is often assumed to be quadratic in 
the total spin operators, and its main term reflects an 
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FIG. 4: (color online) Effective size Dfi^) of the linear su- 
perposition |*}, with Mi = +M and M 2 = -M, for the Fe 8 
spin cluster: the black squares are the computed values, the 
red curve corresponds to the exponential fit. Inset: Same 
function computed for each of the two components |^fc). The 
white squares and the green curve correspond to the com- 
puted values and to the exponential fit, respectively. 



axial symmetry: Hacf — DS^ (with D < 0), where 
z is the principal symmetry axis of the molecule. This 
removes the degeneracy within the ground S multiplet 
of Hexc, giving rise to the characteristic double- well po- 
tential, with the two degenerate ground states M = ±S 
(Fig. Illc)). Hereafter we consider linear superpositions 
of the eigenstates of H = H Bxc + Hacf- These can 
be induced by static terms that don't commute with S z , 
such as the rhombic crystal field Hrcf — E(S% — Sy), 
or a transverse magnetic field. Alternatively, arbitrary 
linear superpositions can be generated by multifrequency 
pulse sequences, that exploit the removed degeneracy be- 
tween the different AM = ±1 transitions [TTj . 



A. The Mni2 molecule 

The Mn 12 molecule represents the first and prototyp- 
ical example of a single- molecule magnet. It essentially 
consists of an external ring, formed by eight Mn 3+ ions 
(each carrying an s = 2 spin), and an internal tetra- 
hedron, formed by four Mn 4+ ions (with s = 3/2). Ex- 
perimental investigation combined with spin Hamiltonian 
calculations allow to group the spins in two sublattices, 
depending on their relative orientation. The sublattices, 
labeled A and B, are formed by the external and internal 
spins, respectively (Fig. [Hb)). The dominant part H Bxc 
of the spin Hamiltonian includes the exchange interac- 
tions between neighboring spins belonging to the same 
sublattice (Ha and H B ) or to different ones (Hab)'- 

H A = JaJZ^-4+i, H B = J B f^sf.sf 

i=l i=\ j=l 

4 4 

Hab = JAB^i-r sf+J^E s ^-( s ? + s ?+i)( 16 ) 

i=l i=l 



where sf = 2, sf = 3/2, s$ = s r , and sf = sf . The 
calculations presented hereafter arc obtained by setting 
[3S]: J A = -64.5 K, J B = 85 K, J AB = 215 K, J' AB = 
85 K. The ground-state multiplet of H Bxc = Ha + H B + 
Has belongs to the S = 10 subspace. 

Before discussing the size of the linear superpositions 
| >F) , we characterize the eigenstates of H Bxc in terms of 
single-spin states and staggered magnetization, in order 
to verify their resemblance to a classical-like spin state. 
The maximally polarized states (Mi = 10) present the 
highest degree of polarization on the external spins sf, 
and more specifically on the even-numbered ones [28] . 
This is reflected by the main populations in the reduced 
density matrices p\ I , which correspond for all spins to 
the states with \rtii\ — Si, with pi(mf = 2) = 0.692, 
Pl (m£ = 2) = 0.845, and Pl (mf = -3/2) = 0.591. The 
main configuration of the M\ — 10 eigenstate is repre- 
sented by the classical-like ground state, with K'I'i = 
10 |vEr?™>*)| = 0.307. The expectation value of the stag- 
gered magnetization and its variance are given by 17.6 
(out of a maximum of J2 X =A b ^2i s f = 22) and 7.00, 
respectively. These values undergo respectively a linear 
decrease and an exponential increase with Mi, that re- 
flects the dependence on the total-spin projection of the 
relevant Clebsch-Gordan coefficients (see Eqs. 12 and 
131. 

The size of the linear superpositions between the 
ground states of H exc is first quantified in terms of the 
quantum Fisher information. The dependence of the size 
Dfi(^) on the total spin projections Mi and M^ of the 
components \^k) is reported in Fig. [2Fa). As expectable, 
the maximum of Dfj(^) is achieved for the maximally 
polarized components (Mi = —A/2 = 10). The size of 
the linear superposition shows a strong dependence on 
I Ml — M 2 \, while it depends only weakly on |Mi + M 2 |. 
At a qualitative level, the same features characterize the 
relative Fisher information, and the corresponding size 
Drfi(^) (Fig. [5Jb)). From a quantitative point of view, 
we note that the maximum value of both the absolute and 
the relative Fisher information is comparable to that ob- 
tained for Fes (see the following Subsection, and Table 
IT]) , in spite of the fact that the sum of the spins gives a 
slightly smaller value in the latter case. 

The size of the linear superposition |>F) can also be 
quantified in terms of the distinguishability by means of 
local measurements of the two components. The proba- 
bility P; of discriminating between |^i) and l 1 !^) with 
single-spin measurements is reported in Fig. [3] for the 
three inequivalent spins, as a function of Mi = — M2 = 
M . On average, the probabilities are slightly lower than 
those obtained for the Fes cluster (see below). In particu- 
lar, we find that only one of the spins (s 2 , green squares) 

reaches values of Pj > 0.99, and only for M = 10. In 
order to estimate the effective size D^m for such case, 
we need to consider different partitions of the spin clus- 
ter, where the spins with least distinguishable states are 
grouped together in larger subsystems A\ , so as to achieve 
values of the probability P^ above the threshold. The 
partition with the largest number of subsystems that ful- 
fils such requirement is defined by the four subsystems 



M = {4 



l' 6 l 



)} (purple squares), with 1 < I < 4, 



and by other four, each formed by an individual spin: 
A l+A = {s&}. As a result, D LM (^,S = 10~ 2 ) = 8 for 
M = 10. 



B. The Feg molecule 

The second most investigated single-molecule magnet 
is commonly indicated as Fes- The magnetic core of Feg 
consists of eight Fe 3+ ions, each carrying an s = 5/2 
spin (Fig. [lja)). Experimental evidence exists that the 
Feg molecule has a ground S — 10 multiplet, resulting 
from a ferrimagnetic ordering of the spins [T| . The spin- 
Hamiltonian calculations allow to further specify the spin 
ordering. This is characterized by six spins (sublattice 
A, blue circles) oriented parallel to each other and to the 
total spin, while the remaining two (sublattice B, red 
circles) are oriented in the opposite direction. The domi- 
nant part of the spin Hamiltonian includes the exchange 
couplings represented in the Fig. flja), and is given by 
H Ex c = H a + H A b, where: 



H A 
Hab 



J a J2 4+i ■ (4-i + 4) + A4 ■ 4, (17) 

4=1 

6 2 

Jab X) X] 4 ■ s f + 

«=5 J=l 

Ji B [4-(4 + 4) + 4 -(4 + 4)] ,(i8) 



with sf = s B = 5/2. A good agreement with the ex- 
perimental data is obtained with the following set of val- 
ues for the exchange constants: J a = 26 K, J' A — 36 K, 
Jab = 200 K, J' AB = 59 K [2S]. Even in the absence of a 
magnetic field, the degeneracy in the ground S multiplet 
is removed by the crystal field. In particular, due to the 
relatively high \D/E\ ratio, the two lowest states in the 
double well are approximately degenerate and with de- 
fined values of the total-spin projection, while a certain 
degree of admixture is induced in the low-|M| states by 
the rhombic term. Hereafter, the two components of the 
linear superposition are identified with the eigenstates of 
Hq = He X c + Hacf (i.e. the part of H that commutes 
with 5*2), that can be coherently coupled by the rhombic 
term Hrcf- 

We start by verifying to which extent the eigenstates of 
Hq approach the classical-like form |\]>™ a:E ) (Eq. [5]). The 
single-spin states corresponding to the maximally polar- 
ized eigenstates M\ = 10 are characterized themselves by 
a strong degree of polarization, where the largest popula- 
tions correspond to the states with the largest spin pro- 
jection: pi(mf = 5/2) = 0.976, pi(m$ = 5/2) = 0.668, 
and pi(mf = —5/2) = 0.643. Remarkable differences 
emerge in this respect between the four external spins 
(x — '- A, 1 < i < 4) and the four central ones. The result- 
ing overlap between the Mi = 10 ground state and the 
main spin configuration is given by |\l/ 1 |\E'5 raa:r )| = 0.687, 
while |(*i|5*|*i) = 28.0. States with a lower M x are 
characterized by less defined single-spin states, and by 
a larger fluctuations of the staggered magnetization S* , 
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FIG. 5: (color online) Probability P ; of discriminating be- 
tween the two components |$i) and ^2) (with M\ = +AI 
and M2 = —M) of the Fe8 spin cluster. The probabilities 
refer to subensembles formed by q spins: red, blue and green 
squares correspond to the spins s$~, sf and sf, respectively 
(q — 1); the purple squares correspond to the four central 
spins S5 , Sq , Si and s 2 (q — 4). Inset: Values of the von 
Neumann entropy S(s*,k) for the three inequivalent spins 
corresponding to each component |$fc). Same convention as 
above for the colors. 



that decreases exponentially with the total-spin projec- 
tion (see the insets of Figs. |i]and[5j respectively). 

We next quantify the size of the linear superpositions 
between states with opposite spin projections: M\ = 
— Mi = M > 0. The two components of such linear su- 
perpositions are increasingly distinguishable in terms of 
S* for increasing values of M, being also (Mi \S* | Mi) — 
(M 2 \S*\M 2 ) oc M. Therefore, the fluctuations of the 
staggered magnetization in |^) grow exponentially with 
M, and so does the effective size Dpi(^) of the linear 
superposition (Fig. H). Such exponential increase e aM , 
with a ~ 0.178, is independent on the specific form of 
the eigenstates, for it reflects the exponential increase 
with M of the Clebsch-Gordan coefficients Aff^ M (see 



Eq. 13). For M — 10, Dfi{^) approaches the maximum 



possible value for this system, ^2 A b Si s f ~ 20- The 
same exponential increase characterizes the dependence 
on M of the relative Fisher information Dhfj(^) (not 
shown), that ranges from 1.148 to 48.707, for 1 < M < 
10. 

We next consider the size of the linear superpositions 
based on the distinguishability of the two components 
l^i) and ^2) by means of optimized local measure- 
ments. The corresponding size of the linear superposi- 
tion is given by the quantity Dlm(^,8) (Eq. 15), here 
with 5 to 10~ 2 . The dependence of Dlm{^>) on M is 
reported in Fig. [5] for each of the three inequivalent 
spins in the cluster. As suggested by the expectation 
values of the above comments on the single-spin states, 
P.' is highest for the external spins belonging to the 
sublattice A (1 < i < 4, green squares), and lies above 
the threshold 1 — 5 for M > 9. The probabilities corre- 
sponding to the central spins (red and blue squares) lie 



instead below the threshold for all the considered cases. 
In order to efficiently discriminate between the two com- 
ponents |\&fc), one needs to perform an ideal measure- 
ment on the whole central core of the cluster, i.e. on 
the spin subensemble A$ formed by sf , sf , Sg , and 
sf (purple squares). Here, the condition Pg > 0.99 
is achieved for M > 9. We note that such condition 
cannot be fulfilled by further partitioning the subensem- 
ble corresponding to the four central spins, for all the 
three inequivalent pairs give P ; < 0.99 (I — 1,2,3) 
for any M (not shown). From the above results, one 
can conclude that D LM (^,S = 10~ 2 ) = 5 for M = 9 
and M = 10, where the optimum partition corresponds 



to A t = {sf} (/ = 1, . . . ,4) and A 5 = {sf 



B „B C .A 



A A\ 
5 1*6 J- 



Therefore, while the larger number of spins doesn't result 
in a larger size Dpi(^) for Mni2 with respect to Fes, it 
does allow to obtain larger linear superpositions in terms 
of £> LM (*,£). 
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FIG. 6: (color online) (a) Effective size of the linear su- 
perposition of the maximally polarized ground states (M = 
±6(sa — Sb)) of the bipartite rings, as a function of sa 
(sb — 1/2). The values of Df/($) (black squares) displays 
a linear dependence on Sb (red line). The maximum val- 
ues corresponding to the four cases are given by the white 
squares, (b) Size of the same linear superpositions quantified 
by Dbfi(^). The red curve represents an exponential fit of 
the computed points. 



III. BIPARTITE SPIN RINGS 

The Mn 12 and Fe§ molecules are relatively complex 
spin clusters, where the polarization of the spins is to 
some extent limited by the competition between differ- 
ent exchange paths. This limits the size of the linear su- 
perpositions that can be possibly generated in these sys- 
tems, both in terms of Fisher information of the staggered 
magnetization and of distinguishability of the two com- 
ponents by means of local measurements. Heterometallic 
spin rings formed by alternate sequences of inequivalent 
spins represent a class of simple spin clusters, that are 
immune to such limitations. The spins sa and Sb < sa 
belonging to the two sublattices, with Na = Ng = N/2, 
are coupled to each other by isotropic exchange interac- 
tion between nearest neighbors: 



N/2 



H+ 



l), 



(19) 



where s^, 2+1 = sf . The antiferromagnetic coupling 
{Jab > 0), combined with the spin difference sa — sb, 
results in the formation an S — N(sa — Sb)/2 ground 
multiplet. A prototypical example of such a spin clus- 
ter is represented by the Mng ring, formed six Mn 2+ ions 
( s a = 5/2) coupled to six organic radicals (sb — 1/2). In 
order to highlight the role played by the inequivalence be- 
tween the two spins, this case will be contrasted with that 
of analogous spin rings, characterized by Na = Nb = 6, 
Sb = 1/2 and sa — 1,3/2,2. We shall focus on the 
most relevant linear superposition, namely that between 
the Mi = +S and M2 = —S eigenstates of the ground 
S- multiplet. 

The dependence of Dpi(^) on the value of sa is 
reported in Fig. [6la). The calculated values (black 
squares) follow a linear dependence (red line) on sa ■ For 
increasing sa, the value of Dpil^) approaches its max- 
imum possible value 6(sa + sb) (white squares). We 
also note that the size of the linear superposition cor- 
responding to the Mng molecule has approximately the 



same value than that obtained for the Fes and Mni2 
spin clusters (Table III) . If normalized to the theoretical 

maximum (J2 X =A b St=i s f = 18), such value is the 
highest amongst the considered molecules. The effec- 
tive size based on the relative Fisher information grows 
exponentially with sa (Fig. [6lb)), and reaches for the 
Mng molecule a value which is much larger that of the 
previously considered molecules (Table IT]). This results 
from the classical-like character of the M\ = 12 states 
(l^rl^™ 01 )] = 0.762), corresponding to small values 
of the fluctuations in staggered magnetization and to 
strongly polarized single-spin states: pi(rn^ = 5/2) = 
0.915 and pi(mf = -1/2) = 0.911. As discussed in the 
following Section, the size of the linear superpositions 
could be further increased if the parallel alignment were 
not only allowed by the absence of antiferromagnetic cou- 
plings within each sublattice, but favored by the presence 
of strong ferromagnetic interactions. 

We finally consider the size of linear superposition 
between the ground states with Mi = — M% = 12 of 
the Mng ring, based on the distinguishability of the 
components by local measurements. This is given by 
Dlm{^,S — 10~ 2 ) = 7. In fact, the two components 
can be efficiently distinguished by single-spin measure- 
ments performed on the spins of sublattice A. Instead, a 
probability above the threshold 1 — 5 cannot be achieved 
with any three spins sf , nor with smaller subensembles. 
Therefore, the partition of the spin cluster into seven sub- 
systems - six including one spin of A each (Ai = {sf}, 
with I — 1,...,6), one including all the spins of B 
(A 7 = {sf , s% , sf , sf , sf, sf }) is the optimum one. 



IV. COMPARISON WITH OTHER SPIN 
CLUSTERS 

The molecular nanomagnets considered so far, Mn.12, 
Fes, and Mng, can be regarded as three representative 
examples of high-spin molecules, with ground states that 
are suitable for generating linear superpositions of rela- 
tively large sizes. It might be useful, however, to high- 
light how such sizes can be further increased, for example 



by engineering the exchange couplings, without modify- 
ing the geometry of the spin cluster. In particular, one 
can consider the case where the dominant part of the spin 
Hamiltonian is still represented by isotropic exchange, 
and the linear superpositions involve the two states with 
Mi = +S and M 2 — —S. Also, the sublattices A and 
B for each of the spin clusters are left unchanged with 
respect to those represented in Fig. [T] We find that the 
states that maximize the fluctuations of the staggered 
magnetization are the ones with maximum values of the 
partial spin sums, and with S^ and S^ antiparallel to 
each other: 

|**=i, 2 ) = \S A =ST X , S B =S% ax ,M=±(S A -S B )), (20) 



where S^% jB 



,Af 



J2i=i s ? 1^9- These states corre- 
spond to the maximally polarized ground states of He X c, 
in the limit of strong ferromagnetic intra-lattice cou- 
plings, and weak antiferromagnetic inter-lattice interac- 
tions: Ja,Jb < 0, Jab > 0, and \Ja\, \Jb\ > Jab- Be- 



sides, 



< &fc=i,2) approach the classical-like states |* 



max \ 
k 



(Eq. p ) for increasing values of the difference Sa — S B ■ In 
Tablelllwe report the sizes of the linear superpositions be- 
tween the eigenstates Ml = S and M 2 = — S for the three 
spin clusters Mni2, Fes, and Mng, along with those of 
the corresponding linear superpositions (|$i) + |$2))/v / 2- 
These latter values (in parentheses) approach the maxi- 
mum possible values that could be obtained within the 
same spin clusters, corresponding to the states \$ ma x)- 

Molecular spin clusters with only ferromagnetic 
interactions can be regarded as the limiting cases of 
the ones considered so far, with Na = N and N B = 0. 
Here, the ground states corresponding to M = ±5, 
with S = Y]a—i Si, correspond exactly to classical-like 
states of the form |\I>™ 0:E ). The sizes of the largest linear 
superposition between M = +S and M = —S are given 
simply by D FI {9) = S and D LM (<5>,5) = N, where the 
latter expression is independent on S. As a remarkable 
example of such a molecular nanomagnet, we mention 
[Et3NH] 2 [Mn(CH 3 CN)4(H 2 O)2][Mn 10 O4(biphen)4Br 12 ], 
Mnio in short [3D]. Such a spin cluster has an S = 23 
ground multiplet, and allows in principle to generate 
larger linear superpositions, with respect to the pre- 
viously considered molecules (Table III). On the other 
hand, we note that the cost of a linear superposition 
(e.g., in terms of the number of microwave pulses) tends 
to increase with the difference between the total spin 
projections, Mi — M 2 . Therefore, molecular spin clusters 
with a ferrimagnetic ground state can in principle allow 
to generate larger linear superpositions with respect to 
ferromagnetic ones, for a given S. For example, the 
values of D FI ($) and D LM (^) obtained for Fe 8 (Table 
IIJ) are larger than the ones that could be obtained within 
a cluster of A = 4 ferromagnetically coupled s = 5/2 
spins: D FI (^) = 10 and D LM (^)A. 

Molecular spin clusters can also exhibit mesoscopic 
quantum coherence in the form of Neel-vector tunnel- 
ing. Prototypical systems where such phenomenon has 
been investigated are even-numbered spin rings, with an- 
tiferromagnetic exchange between nearest neighbors [3"lT - 
[34] . Here, the sublattices A and B have equal sizes, and 





Mni2 Fe 8 


Mn 6 


Mmo 


Dw(») 


14.4 16.5 


16.0 


23.0 




(21.0) (18.3) 


(17.1) 


(23.0) 


£>F/(*fc) 


0.318 0.339 


0.115 


0.0 




(0.140) (0.121) (0.0584) 


(0.0) 


D RFI {^) 


45.4 48.7 


139 


- 




(143) (152) 


(293) 


- 


d lm (\M = io- 2 ) 


8 5 


7 


10 




(10) (8) 


(8) 


(10) 



TABLE I: Size of the linear superposition \i£>) between 
the maximally polarized ground states of the considered 
molecules: Mni2, Fes, Mn6, Mnio. The values in parenthe- 
ses refer to the models of the spin clusters where the two 
macrospins Sa — 2»=i s i an d Sb — ^2iJ[ s f are antiferro- 
magnetically coupled. The values of Dfi{*&) corresponding 
to the states \^ m ax) of Mni2, Fes, Mn6 and Mnio are 22,20, 
18, and 23, respectively. 



are formed by the odd- and even-numbered spins, re- 
spectively. In the simplest case, the easy-axis anisotropy 
tends to enhance in the ground state |^) the weight of 
the classical- like ground states l^^ff 2 )j while quantum 
tunneling between them is induced by an in-plane mag- 
netic field We note however that the application of the 
measures Drfi^) and D F m(^) to such states is not 
straightforward. In fact, the overlap between the actual 
ground state 1^) and \& max ) is only moderate for realis- 
tic values of the physical parameters, and doesn't allow 
to unambiguously decompose |\&) into a linear superposi- 
tion of two components (Eq. IT]). On the other hand, the 
fluctuations of the staggered magnetization alone (and 
thus Dfi{^)) don't allow to distinguish between a meso- 
scopic Schrddinger cat state and one which is quantum- 
mechanical in a broader sense 1191. 



V. CONCLUSIONS 

In conclusion, resorting to different quantum- 
information theoretic measures we have quantitatively 
investigated the size of linear superpositions that can be 
generated in molecular nanomagnets. The prototypical 
single-molecule magnets, namely Mn 12 and Fes, exhibit 
comparable values of the considered measures. Linear 
superpositions of larger sizes can be possibly generated 
in the ring-shaped molecule Mng, at least in terms of the 
relative Fisher information. The absence of competing 
interactions within each of the two sublattices (A and 
B) and a large unbalance between the partial spin sums 
(Sa and S B ) are identified as rule-of-thumb principles 
for maximizing the size of linear superpositions within 
the ground S multiplet and for a given geometry. While 
these criteria are fully met by spin clusters with ferromag- 
netic exchange only, ferrimagnetic systems might allow a 
good compromise between large sizes of the linear super- 
positions and moderate values of the total-spin projection 
difference Mi - M 2 . 
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Appendix A: Maximization of the quantum Fisher 
information 

Hereafter, we maximize the quantum Fisher informa- 
tion J 7 ,/, (A) as a function of the operator A, where X is 
a linear combination of single-spin operators: 



X 



N 

E 



X, 



N 



(Al) 



with n = (n x ,n y ,n z ) and n| = 1. The maximiza- 
tion is performed for the most relevant linear superpo- 
sitions considered above, namely those where |\&i) and 
|^2) have opposite values of the total-spin projection 
(Mi = -M 2 = M > 0). 

We start by considering the case where the versors n^ 
are all identical: X = h ■ J2i=i Sj = n • S. Being the two 
components the time-reversal of one another, the follow- 
ing relations apply: 

(* 1 |A|* 1 ) = -(* 2 |A|* 2 )=7i 2 M,(*|A|*)=0,(A2) 
(*i|A 2 |*i) = (* 2 |X 2 |* 2 ) = (*|A 2 |*)=n2M 2 .(A3) 



In Eq. |A3| we have made use of the assumption that 
Mi — Mi > 2, which allows to neglect the contribution 
to the variance of X coming from the coherence between 
the two components: (5 , i|A 2 |5' 2 ) =0. As a result, the 
Fisher information of the state |\&) reads: 

T i; (X)=2[S(S + 1) - M 2 ]+2n 2 z [3M 2 -S(S +1)}. (A4) 

Here, the equivalence between the x and y directions al- 
lows to set, without loss of generality, n y = and thus 
n x = 1 — n 2 . From the above equation it follows that 
J-^(X) is maximized by the operator 

X = S z forM 2 > S(S +l)/3, 
X = S x for A/ 2 <S(S + l)/3. 



(A5) 



In the case of the previously considered molecules, the 
ground multiplets correspond to S = 10 (Mni 2 and Fes) 
and S = 12 (Mng). These values lead to the conditions 
M > 7 and M > 8, respectively, for having X = S z 
as an optimal total-spin component. For smaller val- 
ues of M, the component with the largest fluctuations 
in 1^) is no longer S z , but S x . We note, however, that 
such fluctuations don't qualify the corresponding linear 
superpositions as good examples of mesoscopic linear su- 
perpositions, for n z — also maximizes the fluctuations 
at the level of the single component. In fact, the vari- 
ance of A corresponding to |^i) is given by the following 
equation: 



which monotonically decreases with n 2 . For linear su- 
perpositions between strongly polarized states (M 2 > 
S(S + l)/3), the operator X = S z thus maximizes both 
the absolute and the relative Fisher information of |\&). 
We next consider a more general class of operators A, 
by admitting different versors fi^, while keeping the con- 
dition Mi - M 2 = 2Mi > 2. The Fisher information of 
the state |\t) reads: 

N 

F^,(X) = 4 ^ {(ni, x nj, x + n hy n hy )(^i\s ltX s J!X \^ 1 } 



ij=l 



N 



(*ikz^|*i)}= ^^(A,*), (A7) 



where we have made use of the relations in Eq. |A2| 
The operator A that maximizes such quantity is de- 
fined by the optimal values of the versor components rii ta 
(a = x, y, z), and depends on the values of the two-spin 
correlation functions that characterize the system eigen- 
states |\&fc). Such problem is in fact equivalent to com- 
puting the ground state of a system formed by N classical 
spins hi, interacting through anisotropic exchange cou- 



plings J^ = {$i\Si 



|*i) and J = (#ik 



l*i) 



As long as the highly polarized components \^&k) are con- 
cerned, there is no competition between the maximiza- 
tion of each contribution Fij(X). In fact, all the versors 
tend to be aligned along z, being 



K*i 



7 l,X°J,X 



*i}| < K*i|*i,»«j,«|*i>| (A8) 



V* 1 (A) = (l-n 2 )[^ + l)-M 2 ]/2, 



(A6) 



for all 1 < i, j < N. Besides, the N spins can be divided 
into two sublattices A and B such that (\&i|si z flj z |\l/i) is 
positive for all pairs of spins belonging to the same sublat- 
tice and negative for all pairs where the two spins belongs 
to different sublattices. The ground state of the classi- 
cal spin system thus simply corresponds to n z ,i = +1 
for all SieA, and n Zt i = —1 for all SitB. This set of 
versors also defines the optimum operator A. We have 
numerically checked that the above conditions for the 
spin-pair correlation functions are fulfilled by the consid- 
ered molecules for large enough values of M. In partic- 
ular, the inequalities coincide with those obtained above 
in the case of identical versors: M > 7 for Fe§ and Mn 12 , 
M > 8 for Mng. We note that the sublattices A and 
B that are defined on the basis of the correlation func- 
tions (^1 1 Si fZ Sj >z |^i) coincide with those introduced on 
the basis of the single-spin expectation values (see Fig. 

As a final remark, we note that any operator A which 
is quadratic in the single-spin operators, such as A = 
J2i j=i(**ij ' s i)(^'ij ' s j); nas identical variance for the 
state |\&) and for the two components |\£fc), provided that 
these are sufficiently polarized (Mi — M 2 > 4, implies in 
fact that the coherence doesn't provide any contribution 
to the expectation value of A 2 ). Therefore, the only 
fluctuations that are captured by J-* (A) are the ones 
that are independent on the linear superposition, and the 
relative Fisher information is always 1. This makes the 
above quadratic operator A unsuitable for quantifying 
the size of the linear superpositions |\&). 
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